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^ , Abstract: We study correlation functions in (0+l)-dimensional maximally supersym- 

^ \ metric U(A^) Yang-Mills theory, which was proposed by Banks et al. as a non-perturbative 

definition of 11-dimensional M-theory in the infinite-momentum frame. We perform first- 
principle calculations using Monte Carlo simulations, and compare the results against the 
predictions obtained previously based on the gauge-gravity correspondence from 10 di- 
mensions. After providing a self-contained review on these predictions, we present clear 
evidence that the predictions in the large- limit actually hold even at small N such as 
N = 2 and 3. The predicted behavior seems to continue to the far infrared regime, which 
goes beyond the naive range of validity of the lOD supergravity analysis. This suggests 
that the correlation functions also contain important information on the M-theory limit. 
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1. Introduction 

It is widely recognized nowadays that gauge theory plays an important role in understand- 
ing quantum gravity and string theory. In this context maximally supersymmetric U(A^) 
Yang-Mills theory in (0-|-l) dimension is of particular interest. This theory was proposed 
as a non-perturbative formulation of M-theory by Banks et al. (BFSS) [1] in an appropriate 
large-TV limit. It is also called, justifiably, "Matrix theory," the term which we will use 
throughout the present paper even for finite N for the sake of simplicity of nomenclature. 
M-theory is a conjectural 11-dimensional theory [2, 3], which was suggested to play a piv- 
otal role in achieving a complete unification of perturbative superstring theories. When the 
11th dimension of M-theory is compactified to a circle or S^/Z2 of radius R, the effective 
10-dimensional theory one obtains at small R is assumed to be perturbative type II A or 
Heterotic Eg x Eg superstring theory, respectively, where the radius is given as R = ggis in 
terms of the string coupling constant gs and the string length ig- These relations, together 
with the string dualities in 10 (or 9) dimensions, connect all the perturbative superstring 
theories. It is therefore naturally expected that deeper principles behind string theory will 
be revealed if one finds a concrete realization of M-theory. 
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The action of Matrix theory coincides with that of the effective super Yang-Mills 
theory describing a collection of N DO-branes [4] in type IIA superstring theory in 10 
dimensions, where the Yang-Mills coupling proportional to the string coupling 

gs . Each DO-brane with the mass 1 / ggis can be regarded as a Kaluza-Klein state of the 
IID graviton carrying one unit of momentum along the compactified 11th direction. The 
proposal advocated in ref. [1] was that Matrix theory in the large- A'^ limit with fixed ^ym 
provides the exact quantum description of M-theory in spatial 9-dimensional transverse 
directions of the infinite momentum frame that is boosted along the circle in the 11th 
dimension. 

Another closely related motivation for this (O-l-l)-dimensional theory is that it can 
alternatively be regarded as a regularized realization [5] of quantum (super) membrane the- 
ory in the light-cone gauge, again in the large-A limit. Supermembranes have also been 
regarded as basic degrees of freedom in M-theory, playing similar roles as strings in lOD 
superstring theory at least in some particular long-distance regime. Matrix theory can 
therefore be a natural basis for approaches to M-theory from the viewpoint of supermem- 
branes, which has been quite an elusive subject for the last few decades. 

Prom these two interpretations, the DO-brancs arc regarded as "partons" for both 
gravitons and membranes in 11 dimensions, and Matrix theory is supposed to describe the 
dynamics of composite systems consisting of a large number of DO-partons. Since the 't 
Hooft coupling constant A = g-y-^N of Matrix theory must be infinitely large in the M- 
theory interpretation, one has to consider the strongly coupled regime of the Yang-Mills 
theory. Understanding the strong coupling dynamics is expected also to shed light on 
important problems in quantum gravity such as a microscopic description of the formation 
and evaporation of black holes in terms of D-branes. In spite of such prospects, it seems 
fair to say that no substantial progress has been made with respect to the dynamics of 
Matrix theory after the end of the previous century. Unfortunately, M-theory itself still 
remains an enigma. 

Although Yang-Mills theory in (0-1-1) dimensions is formally super-renormalizable with 
respect to the UV behavior, there is very little knowledge on its strong coupling dynamics, 
which is essential for understanding the theory in the IR regime. A crucial obstacle is the 
severe IR divergence that comes from loop effects of massless fields in calculating correlation 
functions. While it is expected that some mass scale generated non-perturbatively makes 
the correlation functions finite, perturbation theory does not seem to give us useful insight 
on such a mechanism. In fact, by introducing mass scales through background fields, one 
can perform perturbative calculations, which nicely demonstrate that Matrix theory can 
indeed describe gravitational interactions of DO-branes in agreement with the results from 
supergravity. A remarkable example is the derivation of the general relativistic three- 
body forces [6] from the gauge theory side. This shows that the basic structure of the 
interaction of slowly moving DO-branes is indeed consistent with supergravity in the weakly 
coupled regime. However, for the purpose of studying correlation functions that would carry 
crucial information on the wave functions of DO-brane bound states, we cannot assume any 
preferred background fields. 

On the other hand, the gauge- gravity correspondence has been studied intensively 
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over the past decade as a possible powerful approach to strongly coupled gauge theories 
including QCD and condensed matter systems. The best tested example is the AdS-CFT 
correspondence [7] between 4D A/" = 4 SU(A^) super Yang-Mills theory (SYM) and type IIB 
superstring theory on the special background AdSs x S^. Both sides are characterized by ex- 
act (super) conformal symmetry, which greatly facilitates the analyses. Prom the viewpoint 
of string dynamics, this correspondence was motivated from the open-string/closcd-string 
duality. In particular, if one takes the large- limit with a large but fixed 't Hooft coupling 
A on the SYM side, both the gravitational coupling constant Gio ~ Qg ig and the cur- 
vature (~ {gsN)~^/* £~^) of the background geometry become small on the string theory 
side in the bulk. Therefore, one can study the gauge theory in that limit by using classical 
supergravity theory, which describes the low-energy limit of string theory. 

In view of the general arguments based on the open-string/closed-string duality, there 
seems nothing that restricts the gauge-gravity correspondence to the cases with exact 
conformal symmetry. In particular. Matrix theory is expected to be dual to type IIA 
superstring theory on the background obtained in the near-horizon limit of the DO-brane 
solution [8] . Although neither the action of Matrix theory nor the near- horizon background 
of DO-branes is invariant under conformal transformations, one can transform them back 
to its original form by allowing the coupling constant to transform appropriately. This 
"generalized" conformal symmetry [9, 10] turned out to be useful in formulating the gauge- 
gravity correspondence at the level of operators. Unlike the familiar AdS cases, the dilaton 
and the curvature in the present non-conformal case are not constant in space-time in the 
near-horizon limit of the DO-brane solution. In fact the string coupling becomes stronger 
towards the origin (which corresponds to the IR region in the SYM), while the curvature 
becomes stronger towards the boundary (which corresponds to the UV region in the SYM). 
Thus there is a certain region in which we expect the supergravity description to be valid. 
In order for this region to extend both in the UV and IR directions, N must be sent to 
DO with a fixed but large 't Hooft coupling constant. However, the region of validity of 
this correspondence is slightly different from the BFSS large-A?^ limit that is relevant to 
M-theory. This is not so surprising since the gauge-gravity correspondence is based on 
the DO-brane picture in 10-dimensional space-time, and the space-time dimension must 
somehow be elevated to 11 dimensions for the M-theory interpretation. Let us recall here 
that the dilaton carries information on the effective size of the 11th direction. Therefore, 
one may hope that the predictions from 10-dimensional supergravity contain important 
information on the 11-dimensional theory as one effectively probes the 11th direction by 
approaching the center of the DO-brane solution. This is one of the main issues we would 
like to address in the present work. 

Using the gauge-gravity correspondence, two of the present authors (Y.S. and T.Y.) 
obtained two-point correlation functions in Matrix theory from supergravity on the DO- 
brane background [11]. This work was based on an extension of the general prescription 
proposed by Gubser, Klebanov, Polyakov and independently by Witten (GKPW) [12] for 
the standard AdS case, which enables us to compute correlation functions in gauge theory 
by evaluating the supergravity action as a functional of boundary values of the fields. 
Natural candidates for operators corresponding to the supergravity modes are given by 
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the Matrix theory currents, as suggested from a perturbative analysis with non-trivial DO- 
brane backgrounds [13]. They are single-trace operators with totally symmetric ordering 
inside the trace, which may be viewed as analogs of the BPS-operators in the J\f = 4 SYM. 
The dictionary between these operators and the supergravity spectrum on the DO-brane 
background was given in ref. [11] by matching the scaling dimensions with respect to the 
generalized conformal symmetry. 

For the operators corresponding to the massless modes, which can be described by the 
supergravity approximation, the two-point correlation functions are predicted to take the 
general form 



The supergravity approximation with the near-horizon geometry is valid in the region 



The exponent v in (1.1) turned out to be fractional numbers. They are different from 
those determined by the canonical dimensions of the operators, but are consistent with 
the generalized conformal symmetry. The fact that the exponents are fractional numbers 
suggests that confirming this result from the gauge theory side would require a completely 
non-pcrturbative calculation of the exponents. On dimensional grounds they must be 
independent of the Yang-Mills coupling constant, which has a positive mass dimension. 
It should also be emphasized on the same grounds that the power-law behavior itself is 
quite non-trivial. This suggests that some mechanism for decoupling the non-perturbative 
dynamics from the mass scale of the Yang-Mills coupling constant is at work for this class 
of operators, possibly owing to supersymmetry. 

From the viewpoint of the Matrix-theory conjecture, on the other hand, the power-law 
behavior in the IR regime is an important dynamical signature of the theory [14]. Namely, 
it is related to the existence of a threshold bound state (See ref. [15] and references therein.) 
at each finite N. This is required for consistent identification of DO-branes [4] as partons 
for 11-dimensional gravitons. For unprotected operators, on the other hand, it is natural 
to expect that two-point correlation functions are coupled with the Yang-Mills dynamics of 
mass generation and are given in the form of exponential functions of some dimensionless 
combination of the coupling constant and the coordinate distances. Indeed, such a behavior 
[16] has been predicted for operators corresponding to stringy excited states with large 
angular momentum in the plane-wave limit of the background geometry. 

In the present paper, we study the two-point correlation functions by performing Monte 
Carlo simulation of Matrix theory. The same method [17] has been previously used to 

calculate the internal energy and the Wilson loop in Matrix theory at finite temperature 
by the authors including two of us (M.H. and J.N.) [18, 19, 20]. The results of these works 
are consistent with the predictions from the dual geometry in the bulk. For example, the 
ADM mass of the near extremal DO-brane solution gives the energy density on the gauge 
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theory side as a function of temperature as^ 



IE _ .( \ 14/5 



7-4 TT77 , (1.3) 



assuming that the temperature is identified with the Hawking temperature of the black 
hole. This relation is expected to hold under the condition N^^^l'^^ <^ T\~^/^ <^ 1, which 
coincides naturally with (1.2) under the replacement 1/T — > |t — The Monte Carlo data 
on the gauge theory side not only confirms this prediction but also gives consistently the 
power of the sub-leading term with respect to the inverse YM coupling that are obtained 
from the a'{= £g ) corrections on the gravity (string theory) side [19]. The agreement im- 
plies that the microscopic origin of the black-hole thermodynamics has been accounted for 
in terms of open strings attached to the DO-branes as described by Matrix theory. These 
results encourage us to apply the same method to more general and subtle observables such 
as the correlation functions of local operators. Let us recall that the correlation functions 
are calculated on the gravity side [11] by extracting the response from the DO-brane back- 
ground against perturbations consisting of various different components of gravitational 
waves propagating from the boundary towards the central region. Such observables there- 
fore contain more detailed information on the dynamics of DO-branes than the observables 
that have been studied so far. 

The aim of the present paper is twofold. Firstly, we give an explicit comparison of the 
Monte Carlo data with the predictions for supergravity states [11] and also with those for 
stringy excited states [16]. This provides, for the first time, a highly nontrivial test of the 
gauge-gravity correspondence at the level of local operators in the non-conformal case.^ 
Secondly, we provide evidence for the possibility that the predictions obtained from the 
gauge-gravity correspondence actually hold beyond the naive bounds of the applicability 
of supergravity approximation, both with respect to the size N of the gauge group and to 
the range (1.2). This may be taken as an important piece of information on the behaviors 
of correlation functions in the M-theory regime. We hope that the present work^ would be 
an impetus for revisiting the Matrix theory conjecture from a new perspective. 

The rest of this paper is organized as follows. In section 2, we review the gauge- 
gravity correspondence for Matrix theory and summarize its predictions for correlation 
functions. We hope that this section makes the present paper reasonably self-contained 
and will serve as a useful guide for the readers who are not familiar with the gauge-gravity 
correspondence for DO-branes and Matrix theory. In section 3, we describe our Monte Carlo 
method to calculate the correlation functions and compare the results with the predictions 
from the gauge-gravity correspondence. Section 4 is devoted to a summary and discussions. 



^The power-law behavior (1.3) is also related to the existence of the threshold bound state [21] and 
hence to the power-law behavior of the correlation functions. In contrast, Monte Carlo studies for the 
non-maximally supersymmetric cases [22] confirmed an exponential decrease at low T, which was suggested 
in ref. [21] by assuming the absence of the threshold bound state. 

^See ref. [23] for a numerical analysis based on the discrete light-cone quantization in the (1 -|- 1)- 
dimensional case. 

^Part of the results has been reported briefly in our previous publication [24]. 
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In appendix A, we discuss the physical meaning of the predicted two-point functions for 
stringy operators by considering its spectral representation. In appendix B, we explain the 
details of our Monte Carlo method. In appendix C, we discuss the so-called sign problem 
that appears in Monte Carlo studies of the present system. 

2. Matrix theory and the gauge-gravity correspondence 

In this section we first give a brief introduction to Matrix theory, which was proposed 
as a non-perturbative formulation of M-theory, and discuss the DO-brane solution in 10- 
dimensional supcrgravity, which plays a crucial role in describing the gauge-gravity corre- 
spondence in the present case. Then we review the predictions for correlation functions 
obtained from the gauge-gravity correspondence to the extent that would be sufficient for 
understanding the significance of our main results to be presented in section 3. 

2.1 brief review of Matrix theory 

Matrix theory is the maximally supersymmetric Yang-Mills theory in (0-1-1) dimensions, 
and it can be considered as matrix quantum mechanics. The action is given by 

S = JdtTr (DtXif + [Xi, X,f + i^aDt^Pa - ^ V'a(7^)a/3^^^> V'/?]) , 

9s s gss s ^^^^^ 

where the fields Xj (i = 1, • • • , 9) and (a = 1, • ■ • , 16) are N x N bosonic and fermionic 
Hermitian matrices representing the collective degrees of freedom and their superpartners 
associated with the DO-branes. The theory has a U(A'^) gauge symmetry, and the covariant 
derivative Dt acts as Dt = df — i [A, ■ ] with the gauge field A. The 16 x 16 gamma matrices 
7i satisfy the SO (9) Clifford algebra {7,, 7^} = 2 Sij. The Yang- Mills coupling ^'ym is related 
to the string coupling gg and the string length ig as ^ym — {'^^)~'^9s ^7^- While the action 
(2.1) is written with the Lorentzian signature, we will use the Euclideanized time coordinate 
in actual calculations both on the supergravity side and on the gauge theory side. 

Prom the viewpoint of open string theory, the above action is regarded as the low- 
energy effective theory for the collection of N DO-branes in type IIA superstring theory 
in 10 dimensions [4]. The low-energy approximation is justified in the non-relativistic 
limit, where DtXi and higher derivatives are sufficiently small. The matrices Xi describe 
the lowest modes of open strings connecting the DO-branes, and their diagonal components 
represent the positions of the DO-branes in the 9-dimensional space. The U(l) part, which is 
decoupled from the SU(A'^) part, corresponds to the center of mass motion of the DO-branes. 
The SU(A^) part, on the other hand, represents the relative motions of the DO-branes and 
their interactions through open strings. 

Prom the viewpoint of 11-dimensional M-theory, the DO-brane states represent the 
Kaluza-Klein modes of the graviton supermultiplet corresponding to the compactification 
down to 10 dimensions. The radius of the 11th direction is i? = Qs^s- The BFSS con- 
jecture [1] is a proposal that Matrix theory in the large-A^ limit with fixed Qs describes 
exactly the dynamics of 11-dimensional gravitons in the infinite momentum frame (IMP) 
boosted along the 11th direction with the longitudinal momentum = N/R. The 9 
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spatial directions in the 10-dimensional interpretation now correspond to the transverse 
directions in the IMF. The non-relativistic approximation is justified in this frame since 
P"*" becomes infinitely large in the above limit. Note that the large- A?^ limit considered by 
BFSS is different from the 't Hooft limit, in which g-yu-^ fixed. Let us also note that the 
ordinary weak coupling limit with finite N and small Qs is another way to realize P"*" — > oo 
and hence the IMF. This alternative is usually called the "discrete light-cone quantization" 
(DLCQ). 

The Matrix theory (2.1) is not scale invariant in the usual sense since the Yang-Mills 
coupling constant has a mass dimension. However, it has a generalized conformal symmetry 
(GCS), which can be naturally understood from the viewpoint of the DO-brane dynamics [9]. 
As emphasized in ref. [9] , a deeper motivation for the GCS originates from the space-time 
uncertainty principle in string theory.^ Under the GCS, the coupling constant ^'ym ^ 9s is 
scaled (with the scaling dimension 3) along with the coordinates and the fields as 

t a~H , X ^aX , a^Qs . (2.2) 

Let us recall that the string coupling = e^'t'^ is related to the vacuum expectation value 
{(j)) of the dilaton in the asymptotic region. All the other generators including the special 
conformal transformation can also be defined, and the full conformal algebra is actually 
realized [9]. One can also see that the scaling transformation (2.2) corresponds to the 
Lorentz boost in the IMF since one can combine it with a trivial engineering rescaling 
[t — > a~^t,X — a~^X,is — )■ a~^is) to obtain the boost transformation (t — a~'^t,X — 
X,R~^ a?R) with the IID Planck length ip = g}^^(-s fixed. Thus the GCS is consistent 
with the DLCQ interpretation of Matrix theory. To test the boost symmetry in the BFSS 
interpretation, on the other hand, we need to examine a much more dynamical aspect of 
the theory, since it is related to the scaling behavior in changing N. (See section 2.3.) 
In passing we note that the Hamiltonian of this system is given as 

i7 = -2P-=i?TvQn2-^[X,,X,f + ...) , (2.3) 

where Ilj are the canonical (matrix) momenta corresponding to the matrix coordinates 
Xi. In the M-theory limit, P~^P~ is fixed and P"*" oc is sent to oo. Thus we have to 
consider the spectrum of H which scales as 1/N near the threshold ~ 0. In terms of 
two-point correlation functions, this amounts to considering the far IR region in which the 
time coordinate scales linearly in N. 

2.2 DO-brane solution and its near-horizon limit 

The crucial point for the gauge-gravity correspondence is that the DO-brancs, which are 
described by the action (2.1) from the viewpoint of open strings, can also be described as 
a classical background in 10-dimensional supergravity from the viewpoint of closed strings. 

*For a comprehensive discussion on the idea of the space-time uncertainty relation in string theory, see 
ref. [25], which also contains an extensive list of references. 
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The solution [26] corresponding to a stack of A'' DO-branes is given in the string frame as 

ds'^ = -h-^/^dt^ + h}''^(dr'^ + r'^d^lfj , 

e^ = gsh^'\ Ao = ~(h-^-l), 

^=1 + 4, q = mT,\gsN)tl , (2.4) 



where ^ and j4o represent the dilaton and the RR gauge field, respectively. String theory 
(or supergravity) on this background in the near horizon limit ^ — > ^ is conjectured to 
be equivalent to Matrix theory [8, 9] under certain conditions explained below. Unlike the 
case of AdSs x S^, which corresponds to the near-horizon limit of the D3-brane solution, 
both the dilaton and the curvature depend on the radial coordinate. From the explicit 
form of the background, one can obtain the range of validity 

(5.iV)^/3iV-^/2i « L « (5^iv)V3 ^ (2.5) 

in which the weakly coupled supergravity description is reliable. The first inequality is 
necessary for the string coupling to be weak e"^ ^ 1, while the second one is necessary for the 
curvature to be small in the string unit. Taking the near-horizon condition j- <^ (ggN)^^'^ 
also into account, we find that a wide range of validity can be obtained if we consider the 
't Hooft limit — >^ 00 in the strongly coupled region, namely gsN 2> 1. 

Let us note here that the near-horizon limit of the DO-brane metric can be obtained 
by a Weyl transformation from the AdS2 x metric as 



ds^ = g2/7/,3/i4 W f +(de^ + cos2 edi;^ + sin2 edQl) \ . (2.6) 




Here the Poincare coordinate z is related to r by 

,= lqi/\-5/2 ^ (2.7) 
o 

and we have chosen a special parametrization of S®, which will be useful in section 2.4. 
Given this metric (2.6), the boundary (r 00) is formally located at z ~ 0. One should 
keep in mind, however, that the near-horizon condition <C g in the original metric (2.4) 
puts a natural cutoff around z ~ g^/^. 

As in the standard AdS-CFT correspondence, the time scale At, with which we probe 
the gauge theory, corresponds to the radial scale Az ~ At from the bulk viewpoint, since 
the time is common to both the bulk and boundary theories. (We will see this more 
explicitly in the following calculations.) This relation enables us to convert the region of 
validity (2.5) of the supergravity analysis into 

A-i/-^ « At « A-i/3^i°/2i (2.8) 

in terms of the time separation in the (0 -|- l)-dimensional gauge theory, where we have 
used 

A = g^^N = {27r)-^gsm-^ . (2.9) 
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This region (2.8) with milder N dependence A^^^/^^ does not overlap on the IR side with the 
M-theory regime At ~ A~^/^A^ mentioned below eq. (2.3). However, the Monte Carlo data 
presented in section 3 suggest that the supergravity results may be valid even in the IR 
region At > A-V37V10/21 beyond the naive region of validity (2.8) at least for supergravity 
modes. 

We emphasize that the whole near-horizon background including the Wcyl factor, the 
dilaton and the RR 1-form A^dt (up to an irrelevant constant term for the 1-form) is 
invariant under the same generalized scaling symmetry that appears on the gauge theory 
side, where we transform the parameter gg o?gs together with the coordinates t a''H, 
z — >■ oT^z. We can therefore classify the linearized waves in supergravity around the near- 
horizon background according to the transformation properties under the GCS. Although 
the behavior of massive fields such as stringy excited modes is qualitatively different from 
the massless case, we can still obtain constraints from this scaling property. See also ref. [27] 
for a discussion on the extension of this symmetry to the stringy level in a slightly different 
but related context. 

2.3 predictions for supergravity modes 

Let us now review the calculation of correlation functions from the supergravity side [11]. 
The basic idea is to apply the general prescription [12] for the gauge-gravity correspondence, 
which has been tested by many nontrivial examples in the standard AdS case. The pre- 
scription, formulated in the Euclideanized space-time, states that the correlation functions 
of Matrix theory operators 0^{^) can be calculated by evaluating the bulk supergravity 
action Ssg as a functional of the values of the bulk fields at the boundary as 



^ = Li:i!dthi{t)oi{t)\ ^2.10) 



The supergravity fields classified by an appropriate basis {/} are required to take a 
fixed value = hQ{t) at a regulated boundary z = e. In our case it is natural to assume 
that the boundary is located at the end of the near-horizon region e = g^/^ (or r = 
^i/7-j_ r^Yic fields are assumed to vanish in the central region {z — )■ oo). This condition^ 
enables us to circumvent the singularity at the center, and plays a crucial role in the 
following calculations. The right-hand side of (2.10) represents the generating functional 
for connected correlation functions of Matrix theory operators 0^{t), each of which couples 
to hl{t). 

In ref. [11], the complete spectrum of the linearized fluctuations around the near- 
horizon DO-brane background was obtained. After diagonalizing the fluctuations and de- 
composing them into S0(9) spherical harmonics, each eigen-mode of the bosonic fluctua- 
tions is described by the (Euclidean) effective action 

<9efr = ^ J dt dz z ({dohf + {d,hf + ^/i^ j _ (2.II) 



^From a mathematical point of view, it is needed for the kinetic operators of hnear perturbations around 
the singular background to be self-adjoint near the center. 
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The field h is normalized in such a way that they are dimensionless and the action has 
an overall factor of where ^ £g represents the gravitational coupling. The 

constant u is given by 

9 7 

(2.12) 



2 . 7 
I' = - I + -n 

5 5 



where i is an integer which represents the SO (9) total angular momentum, and n is an 
integer which depends on the type of the field h {—1 < n < 3 for bosonic modes). In 
Table 1 we list the modes which we study by Monte Carlo simulation in section 3 (They 
are denoted in ref. [11] as V2, 02 and S3, respectively.), and describe their properties such as 
the transformation property on S^, the value of u, and the range of £. The origin of these 
fields arc cither the graviton qmn or the three- form potential Cmnp in IID supergravity. 
The corresponding operators in Matrix theory will be specified later as (2.18)-(2.20). For 
a complete list of the spectrum, see ref. [11]. 



mode 


IID origin 


g8 repi-esentation 


V 


range of I 


operator 




Cmnp 


anti-sym. 2-form 




i > 1 


Jl 


a2 


Gmn 


vector 


¥ 


i>2 




S3 


Gmn 


scalar 


i(2£-7) 


i>2 


rp+ + 



Table 1: Supergravity modes and the corresponding operators in Matrix theory. 

Let us compute the two-point correlation functions corresponding to the supergravity 
modes following the prescription in ref. [12] (in particular, the first paper). Imposing the 
boundary condition /io(i) = / dpe'^^^ h{p) at z = e, one can write the solution to the wave 
equation obtained from the effective action (2.11) as 



h{t, z) 



f 



(2.13) 



where K^{pz) represents the modified Besscl function of the second kind, which decays 
exponentially as z — )• 00. Upon substituting the solution (2.13), the action (2.11) can be 
reduced to a boundary term. Setting the boundary at 2; = q^^"^ , and expanding the action 
in powers of pq^l'^ 1 we obtain 



dt 



zhdzh 



z=oo 



= J dph{p) h{-p) (analytic in p) + c {p q^/''^ 1 1 + O (^(p q^^'^f^ | 



(2.14) 



where c is a non-zero constant factor.^ The terms analytic in p correspond to local di- 
vergences (delta functions and its derivatives) in the coordinate space. We ignore these 
terms since they do not contribute to the correlation functions at finite separations, and 



®This factor gives the normalization of the two-point function, but we leave it as a free parameter in the 
following analysis. The correct normalization of some operators, such as the energy-momentum tensor, can 
be fixed by computing two- and three-point functions and requiring the consistency with the Ward identities. 
This was done in ref. [28] using an expansion in the Fefferman-Graham form. We thank K. Skenderis for 
discussions on this point. 
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furthermore we keep only the leading terms which are most relevant to the IR property of 
the correlation functions. From the relation (2.10), we obtain the connected part of the 
two-point function at large separations as 

0(tmt')) = l.pl,,U"'e-^"' ' s (2.15) 



h{p) h{jp') 

/,,e«'-"),-.J,^A^, (2.16) 



which leads to the result (1.1). 

In the last step of (2.16), we have to recall that in order to have a well-defined Fourier 
transform J^^dp F{p) e^^^ , the integrability condition J^^dp F{p) < oo must be 
satisfied. As usual, the nonintegrability for \p\ oo can be treated as distributions, which 
have appreciable supports only at short distances with respect to time separation and hence 
can be ignored for studying the IR behaviors. On the other hand, since our correlation 
functions behave as ~ \p\'^'^ at small p, the inverse Fourier integral for the operators with 
V < — I {e.g., defined below) is divergent at p ~ 0. In these cases, we would have to 
invoke analytic continuation with respect to the exponent v in defining the Fourier integral. 
The resulting formula is well known (t / 0) , 

]= dp \p\^^ = 22-+V2 + . (2.17) 



Since the predictions are not trustable in the UV region, we should use this formula only 

for sufficiently large 

In the above calculation the information near the boundary plays a decisive role in 
discriminating various operators since the boundary condition at the center is chosen such 
that the fields should vanish. All we need to know is the way in which the wave function 
is "reflected" from the central region, where the wave function (2.13) decays exponentially 
as z — > oo. This vanishing behavior near the central strong coupling regime may be the 
reason why our result (2.16) can be valid in the IR region beyond the naive region of validity 
(2.5). On the other hand, it should also be emphasized that while these calculations are 
performed totally within the lOD picture, the waves reflected back to the boundary can 
elicit important information by sinking into the near-central region during the reflection, 
where the dilaton grows and hence the 11th direction begins to open up indirectly. The 
further we go to the IR region of the gauge theory, the longer the corresponding waves stay 
near the central region. Consequently the reflected waves in the IR region at the boundary 
can in principle store richer traces of the 11th direction than in the UV region. 

Let us now describe how to flnd the operators in Matrix theory which correspond to 
the supergravity modes discussed above. For that purpose it turns out to be useful to 
study the currents and their moments in Matrix theory within perturbation theory around 
appropriate backgrounds for the matrix variables following ref. [13]. These operators are 
written in a single trace form, and are coupled to gravitational perturbations. One can 
identify them by analyzing the one-loop effective potential between diagonal blocks, which 
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represents gravitational interactions. The moments are constructed from the basic opera- 
tors by inserting the Xi fields inside the trace and symmetrizing their ordering. Following 
the IID light-cone notation in ref. [13], let us define the operators^ 

= 'TT Cij;h-ie Tr (FijXi^ ■ ■ ■ Xj^] , (2-18) 

= -V Ci;h-ie ^ ((A^O^n • • • Xi,) , (2.19) 

ys^s ^ ' 

where Fij = —i [Xi,Xj\/lg ^ which we will identify with the supergravity modes V2, a^^ 
in Table 1, respectively. We have defined the dimensionless matrix variable Xi = Xijc^l'^ ^ 
which naturally corresponds to the dimensionless combination jq appearing in the DO- 
background in the bulk theory. We have also assumed the same global prefactor 1 / ggig as 
in the DO-brane action (2.1) so that the engineering dimensions of these operators are now 
1. It turns out that this normalization is necessary for matching the bulk modes and the 
gauge-theory operators with respect to the GCS. 

In order for the operators to have definite SO (9) angular momenta, the constant co- 
efficients C's must satisfy the following conditions. The coefficient Ci^,,,i, in (2.20) should 
be totally symmetric, and it should also be traceless under contraction of any two indices. 
The coefficient Ci-i^,„i, in (2.19) should be totally symmetric, and it should also be traceless 
with respect to the indices ii . . , and anti-symmetric under the exchange of i and any 
of The coefficient Cij-i^,,,ii, in (2.18) is totally symmetric, and it should also be 

traceless with respect to the indices ii. . .i£, and anti-symmetric in i,j as well as for the 
exchange of i or j with any of ii . . . . 

The guiding principle in relating the gauge-theory operators with the supergravity 
modes is that the generalized conformal dimensions should match between the bulk fields 
and the gauge-theory operators. Let us define the generalized conformal dimension A of a 
gauge-theory operator 0{t) such as (2.18)-(2.20) by the scaling property 0{t) — t- 0'{t') = 
p^O{t) under t ^ t' = p~^t, gs ^ 9s = 9s- Having in mind the calculation from the 
gravity side based on the effective action (2.11), we assume that the two-point functions 
of these operators obey the power-law behavior, and that the only length scale allowed is 
g^/^ apart from the gravitational constant Qg ig appearing as the overall coefficient. Then 
the GCS along with the usual dimensional analysis fixes the behavior of the correlation 
function with the above normalization to have the general form 

(0{t) 0(0)) ~ |t|-(7A+12)/5 _ (2.21) 
9s 

Comparing this with eq. (2.16), we get A = — 1 -|- ^v. We can see from Table 1 and 

(2.18)-(2.20) that V2-, 02, S3 correspond to , , T^^ , respectively, since the generalized 

'^The operators given in ref.[13] are of the form J^;f^, j^ = ^j^Tr {FijXi^ ■ ■ ■ Xi^). Here we have taken 
special combinations j/, which transform irreducibly under SO (9). The same remark apphes to the oper- 
ators and T/'*' as well. 
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conformal dimensions A coincide if we assign the canonical value 1 to Xi as suggested from 
eq. (2.2) without any anomalous dimension. (See ref. [11] for a complete dictionary between 
the supergravity modes and the Matrix theory currents.) Note that each of Xi = Xi/q^f^ 
contributes 1 — | = | to the scaling dimension, and hence | to the index v in accord 
with (2.12). It is remarkable that the scaling properties can be explained with such simple 
assignment of dimensions with respect to the GCS. We emphasize that from a purely gauge- 
theoretical point of view, the appearance of the factor g^/^ is genuinely a dynamical effect, 
which is difficult to understand without invoking the dual gravity theory. 

As has been mentioned below eq. (2.2), the generalized scaling transformation used 
above is essentially equivalent to the boost transformation along the 11th direction, where 
the longitudinal momentum = N/R is scaled by treating 1/R as a variable with fixed 
A^. In the M-thcory interpretation of the gauge theory, on the other hand, we fix R or gg 
instead and increase N to realize the IMF. It is therefore interesting to examine the above 
general form from this point of view. We consider the transformation N ^ pN together 
with t ^ pt. Then the two-point functions (2.16) scale as 

{0it)0{t')) ^ -^ ^^^^"['^^^^^ , (2.22) 

from which we obtain the weight c^m = — 6z^/7 for each operator under the M-theory boost. 
In terms of the IID light-like coordinates, the exponent v is expressed as 

v=l{l-n^ + n-) + li , (2.23) 

5 5 

where n+ (n_) is the number of upper (— ) light-cone indices in the operator. Thus the 
weight dyi is given by 

dM={l + ^) {n+ - n_ - 1) - + ^) ^ . (2.24) 

This should be compared with the kinematical weight for the boost c^j^j™'* = (?t-+ — n_ — 
1) — where the term —1 in the parenthesis comes from the fact that the currents are 
supposed to be integrated over the x~ direction, and the factor of 1/7 in front of i comes 
from our normalization of transverse fields Xi = Xj/g^/^. The weight cIm found from the 
gauge-gravity correspondence is indeed determined solely from the 11-dimensional index 
structure of the operator, but we observe interesting anomalous factors. It is therefore 
important to clarify whether the behavior (2.16) continues to be valid in the M-theory 
regime corresponding to the far IR region |t — t'| oc A^. Our Monte Carlo data presented 
in section 3.2 seem to suggest that it does. Then, it would be interesting to clarify the 
meaning of the anomalous behavior indicating that the transverse size is compressed by 
the factor of p^/^ under the boost. See ref. [14] for further considerations on this issue. 

2.4 predictions for stringy excited modes 

In this subsection we extend the calculation of correlation functions to operators that 
correspond to stringy excited modes on the gravity side. We note first that this is a highly 
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nontrivial issue since one has to somehow generahze the GKPW prescription, which is based 
on the supergravity approximation. In the most general case, one would have recourse 
to superstring field theory for closed strings, which, however, has never been formulated 
successfully in non-trivial classical backgrounds even at the linearized level. This forces 
us to formulate the gauge-gravity correspondence in the first-quantized picture. Here we 
review a bulk analysis including stringy excited states [16]. The analysis treats a single 
string moving in the direction with J units of angular momentum with a small but finite 
excitation number. We can then use a semi-classical approximation for the center-of-mass 
motion of the string by taking the large- J limit, which is often called the plane- wave limit. 

Following the well-known work in the AdS5 x case [29] , we identify the string states 
in the plane-wave limit with the gauge-theory operators which are the counterparts of the 
so-called BMN operators in the H = 4, SYM. As a typical operator, let us consider 

J 

Oij^^ = ^e^^'^''I^Tv(^XiZ^XjZ^-^^ , (2.25) 

fc=0 

with Z = + iXg, where the X^ and Xg are taken to be the directions of classical 
trajectories in S^. The "impurities" Xi and Xj are the fields in the other transverse 
directions (i, j = 1, • • • , 7). The operator (2.25) corresponds to a state with two oscillators 
at the n-th level being excited. Supergravity modes discussed in the previous subsection 
correspond to the operators with n = 0. We can also consider excitation of oscillators in 
the t direction, which corresponds to an operator with the covariant derivative Dt inserted 
as an impurity. 

The basis of this analysis is the prescription [30] introduced originally to solve some 
puzzles in the original BMN proposal. More recently, this method^ has been applied sys- 
tematically to the evaluation of various correlation functions including three-point functions 
in the AdSs x case [32, 33]. Instead of computing the energy of strings in the global AdS 
and identifying them with the scaling dimensions as in ref. [29], one computes the transi- 
tion amplitude along a trajectory in the Euclidean AdS (with a non-trivial Weyl factor in 
the DO case) connecting two points at the boundary. Note that the real trajectories in the 
Minkowskian AdS, in general, do not reach the boundary, while the tunneling trajectories 
in the Euclidean AdS start and end at the boundary. The amplitudes along the tunneling 
trajectories can then be identified with the gauge-theory correlation functions. This gives 
us a definite way to incorporate stringy excitations. In the present DO case, in particular, 
it is crucial to use this approach since the original BMN approach would require us to 
consider the strings moving in a singular region near the center. That the tunneling tra- 
jectories do not reach the singularity at the center corresponds to the boundary condition 
in the supergravity analysis that the wave function should vanish at the center. 

We now review the key steps of this approach. For simplicity, we set a'{= = 1 
throughout this subsection. Although the original analysis was given for general Dp-branes 
with p < 5, we restrict ourselves to the p = Q case, which is relevant to the present work. 
To explain the idea in a simple setting, let us first start by considering a geodesic for an 

®For a compact review of this approach, see ref. [31]. 
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ordinary massive particle on the Euclidean AdS2, which can be obtained by minimizing 
the action 

S = m dTy/g^udrXf'drX'' , (2.26) 

where g^y is the metric for Euclidean AdS2 given by ds^ = {df!^ + dz'^)/z'^ . If we set the 
mass^ to m = |j, we can regard this action (2.26) as that of a Kaluza-Klein particle 
with J units of angular momentum on the compactified with the time coordinate being 
Euclideanized. The geodesic of our interest is a half circle^^ 

i 

t = £tanhT, z= — - — , (2.27) 
coshr 

which connects two points on the boundary separated hy \tf ~ ti\ = 2L One finds from 
(2.27) that the radial cutoff {z > e) is related to the proper time cutoff (|t| < T) as 

e- 2le-^ = |i/ -ii|e-^ . (2.28) 

Substituting this solution into the action (2.26), we get = m J^rpdT = 2mT. The 
amplitude in this geodesic approximation is then proportional to 

g-5ei ^ g-2mT ^ f ' . (2.29) 



This gives the large- J limit of the two-point function for the operator Tr(^Z'^^, which is 
included in (2.20) as a particular case. The J dependence is indeed consistent^^ with the 
result (2.16) obtained using the wave picture for the operator T^'^ with i = J. 

In order to treat more general operators, we need to embed the above point-like classical 
trajectory in string theory. For that purpose we first perform a double Wick rotation [30] 
t —7- —it, tp —7- —iip in the DO-brane metric (2.6). One can check that the classical trajectory 
discussed above for the Euclideanized AdS2 case agrees^^, in the two-dimensional subspace 
of time and radial coordinates, with a point-like ground-state solution for the standard (so- 
called "Polyakov-type" ) worldsheet action for strings around the DO-background satisfying 
the Virasoro conditions. The solution is independent of the worldsheet coordinate a and 
is given by eq. (2.27) and ijj = |r. The period of a is chosen to be < a < 27ra with a = 
5 J/ (2q^/'^) so that the canonical momentum for -0 is equal to J, i.e., J = -^(f'^'^ J daip. The 
semi-classical amplitude (2.29) can then be reproduced by evaluating the Routh function 
for the worldsheet string action along this solution. 



^The factor of 2/5 in the mass m is due to the fact that in the DO-background, the radius of is 5/2 
in units of the AdS radius. Note also that we have ignored the J-independent term in the mass m since we 
are dealing with the large-,/ limit. 

^'^Here we take the coordinates {t,z) to be dimensionless. In order to retrieve the original coordinates, 
which scale like (2.2) under generalized conformal transformation, one has to multiply them by q^^"^ . 

^^The J-independent term in the exponent can be obtained by computing the zero-point energy of the 
superstring as described in the second paper of ref. [16] , and it also agrees with the result from the GKPW 
prescription. 

^•^The difference of the Weyl factor in the background metric can be absorbed, along the classical trajec- 
tory, by a time-dependent reparametrization of t. 
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The next task is to expand the worldshect action around this classical solution to the 
quadratic order in fluctuations. The quadratic terms are given as^^ 



S^'^ = ^ldrl^^"da [x' + t-\t) x'' + mlx'' + yf + r-^(r) + yf) , (2.30) 
where r(r) represents the trajectory in the radial direction 

which is obtained from (2.27) using (2.7). The fields x and yj (? = 1, • • • , 7) in (2.30) 

represent fluctuations within the {t, z) direction and along the direction, respectively. 
These fields have mass rrix = 1 and niy = 2/5 due to the curvature around the geodesic. 

The supergravity modes can be obtained by restricting ourselves to point-like con- 
figurations, which corresponds to setting x' = = in the worldsheet action (2.30). 
Then the problem reduces to that of ordinary harmonic oscillators with the Hamiltonian 
H = m^a'^a+ ^ written schematically for each of the x and yi oscillators. The amplitude 
from r = — T to r = T is given by 

/ \ 2TO(ata+l/2) 

^_2HT ^ g-2m(ata+l/2)r ^ / f \ _ ^2.32) 

An insertion of (i = 1, • • • , 7) into the ground state operator Tr^Z"^^ corresponds to 
exciting one of the yi oscillators with m = m^, which increases the power oi\tf —ti\ in the 
two-point function by 2my = 4/5. Similarly, an insertion of Df corresponds to exciting the x 
oscillator with m = rux, which increases the power by 2mx = 2. These rules are consistent 
with the results from the wave-function analysis shown in Table 1. For instance, let us 
consider the operator T^~^ with £ = J and with all of the transverse directions confined 
in the Z-plane. Inserting Dt to it creates an operator T^^ with £ = J — 1. The value of 
ly increases by 1 as one can see from Table 1, which implies that the power increases by 
2 in accord with the conclusion obtained from (2.32). Note that, in the present tunneling 
picture, the time direction t of gauge theory is transverse to the r-direction near the 
boundary. As a consequence, we can naturally identify Dt with oscillation in a transverse 
direction, which supports the consistency of the present analysis. In the Minkowskian 
picture, this correspondence is unclear. 

In order to study stringy excited modes, we need to allow x' and y'^ to be non-zero 
in eq. (2.30). Since the x''^ term and the y'-^ term have r-dependent coefficients, we have 
to solve the problem of quantizing oscillators with time-dependent masses. Let us explain 
how this can be done by taking one of the yi oscillators as an example. 



^^We have fixed the worldsheet metric as y/hh'^'^ = {y/hh'^'')~^ = r^^'^{T), as in the second paper of 
ref. [16]. With this choice, m^, my are constant, but x'^, y'^ get r dependent coefficients. If we use the 
conformal gauge as in the first paper of ref. [16], 771^, my become r-dependent. Final results do not depend 
on the gauge choice. 

^''impatient readers may skip these details and go directly to the final result (2.43), which is valid for 
states with high wave numbers. 
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The equation of motion for the n-th Fourier mode in a is given by 

^ y{r) = rn^ir) y{r) , m^r) = r-\r) J + | • (2.33) 

Let the functions f±{T) be the solutions to cq. (2.33) with the boundary condition /±(r) 
as r — 7- iboo, respectively, which arc normalized by /+/- — /-/+ = 1- Let us regard ^(t) 
now as an operator and denote its conjugate canonical momentum by p{t) = i dryir) 
satisfying the canonical commutation relation [y,p] = i. Then we can write the solution to 
the time-evolution equation (2.33) as 

yiT) = U{T)a + f.iT)a^ , (2.34) 

where a and are r-independent operators satisfying [a, a}] = 1. Due to the time reflec- 
tion symmetry of the problem, we have y^(r) = y{—T), which corresponds to the reality 
condition in the real-time formulation, and hence f±{T) = fzf{—T) since f±{T) are real in 
our case. 

Using the Hamiltonian H{t) = ^ (^p^ + m'^{T)y'^^ , we can define the transition ampli- 
tude between an initial state at r = —T and a final state at r = +T along the tunneling tra- 
jectory. The corresponding "Euclidean" S-matrix is given by S{T) = Texp — J^j. dTH{T 

Since the Hamiltonian is quadratic in a and , the S-matrix takes the form 

SiT) = N{T) : exp Q ^(T) (a^f + B(r) a^a + \ A{T) j : , (2.35) 

where the coefficients are expressed in terms of f±{T) as 

A = -l(l±^ + i±^] , N' = i + B = -. . (2.36) 

We can rewrite S{T) as 

S{T) = N{T) exp (^-O , (2.37) 

where b, satisfying [6, b^ = 1 are related to a, through a T-dependent Bogoliubov 
transformation, and is given by 

coshn = -i^l + B + ^-^j . (2.38) 

In fact, with our boundary conditions for /±(r), it turns out that ^4 ~ (1 -|- S) ~ 0, 6 ~ a 
and J\f ^ {1 + By/"^ when T — >■ oo. Thus the final result takes the simple form 

t , 1 /o / ■ \ -{a^a+1/2) 

S{T) ~ (1 + Bf^+^/^ = (2f_{T)f_{T)j . (2.39) 

Let us then discuss the behavior of stringy excited modes. Near the boundary (r — >■ oo), 
the constant term ^ in cq. (2.33) is dominant over the time-dependent string mass term 
(r^'^(r)n^/a^). This results from the fact that the curvature of the background in the 
string unit becomes strong as r — )■ oo. However, when we study the large distance behavior 
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\tf — ti\ — > oo, the geodesic passes deeply through the interior region, where the string 
mass term dominates. In this situation it is difficult to determine the functions /±(t) for 
general excited states. We therefore restrict ourselves to excited states with higher wave 
numbers satisfying n/a ~ nq^l'^ /J S> 1. This allows us to neglect the terms other than 
the string mass term, considering that the regulated boundary is placed at the end of the 
near horizon region (r(r) ~ 1). Therefore the equation of motion can be approximated as 

y = ^y, (2.40) 



where r' is defined by (It' /dr = r~^/^, and y = r^/'^y. The solutions to this equation are 



In terms of r', the boundary (r oo) is reached at finite time Tb, which is evaluated as 

^ I r ^ p/5 . (2.42) 



Plugging this into eq. (2.39) with T = Tb, we obtain the S-matrix as 
5(T) ~exp|-^rb(^ata + ^ 

~ exp |-c(n, J) (^a^a + - t^r^/^j ^ (2.43) 

where c{n, J) is proportional to \n\/J. In the last equality, we have substituted i in Tb by 
q~^/'^\tf — ti\, where U and tf represent the time coordinates on the boundary at both ends 
of the trajectory. Note that the combination q^^^\tf — tip^^ in the exponent is invariant 
under the generalized conformal transformation. In the case of general Dp-branes (p < 5), 
this combination takes the form q^/^^~P^\tf — ti\'^^~P^/^^~P\ 

Note that, in the conformal case p = 3, the power of\tf — ti\ in the exponent vanishes. 
Therefore, the power-law behavior of correlation functions is modified for stringy operators 
[32] only by the emergence of anomalous conformal dimensions, which depend both on the 
Yang-Mills coupling constant and on the wave number n. In contrast to this, for the present 
non-conformal case, we found the peculiar exponential behavior (2.43) in the IR region. 
This is an interesting prediction, which is worth being tested. We consider that the simple 
form \n\/ J of the coefficient in front of the GCS invariant combination q^^^\tf — tjp/^ on 
the exponent may be valid only in the present limit n/a S> 1. In general, it could happen 
that the GCS invariant combination ^q^^^\tf — U]^^^^ appears on the exponent with some 
power rj ^ 1. In section 3.3 we will provide some evidence that the form (2.43) is actually 
valid for stringy operators even with small J and n. Since the predicted behavior is quite 
different from the standard behavior of a massive theory, we discuss the meaning of this 
form based on its spectral representation in appendix A. 
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3. Monte Carlo calculations on the gauge theory side 

In the previous section we have reviewed the calculation of two-point correlation functions 
in Matrix theory from the gravity side based on the gauge-gravity correspondence. In 
particular, the power-law behavior (2.16) was predicted for the operators (2.18)-(2.20), 
which correspond to the supergravity modes as summarized in Table 1. We also have a 
prediction (2.43) for the type of operators (2.25) corresponding to stringy excited modes 
in the large-J and large-n limits. In this section we calculate these two-point correlation 
functions directly on the gauge theory side by a Monte Carlo method. We will see that our 
results for A'" = 2 and N = 3 already show striking agreement with the predictions from 
the bulk side. 

3.1 putting Matrix theory on a computer 

We consider the Euclidean version of the action (2.1), which is given by 

^ = ^ / dtTr Q {DtXif - 1 [Xi,Xj]' + \i^aDt^l^a -\^alf[Xi,i^^]) ■ (3.1) 

Since the coupling constant Qym '^^^ absorbed by appropriate rescaling of the fields and 
the coordinate t, we set the 't Hooft coupling A = g-yu^ ^° unity without loss of generality. 
Then the strong coupling limit amounts to the IR limit. Since the U(l) sector of the 
U(A'') theory is decoupled from the rest, we actually study the SU(A/') theory. 

In order to put the theory on a computer, we need to make the field degrees of freedom 
finite by introducing UV and IR cutoffs in the i-direction. The IR cutoff is introduced by 
compactifying the t-direction to a circle of circumference /?. Since we are not interested in 
the thermal properties of the system in this work, we impose periodic boundary conditions 
on both bosons and fermions so that the supersymmetry is not broken by finite /? effects. 
The UV cutoff is introduced by a sharp cutoff in the momentum space following the pre- 
vious works [17, 18, 19, 20, 24]. This has the following advantages over the conventional 
lattice approach [34]. Firstly the breaking of supersymmetry is milder than in the lattice 
regularization, and it restores quite fast as the UV cutoff is removed [17]. Secondly one can 
implement the Fourier acceleration [35] without extra cost, which makes the simulation 
much more efficient. 

When we introduce the momentum cutoff, we have to take care of the gauge symmetry. 
In the case of one-dimensional theory in general, one can fix the gauge completely {i.e., 
choose a unique representative for each gauge orbit non-perturbatively) in such a way that 
higher momentum modes are more suppressed by the kinetic term in the action. With this 
gauge choice, one can also show by simple power counting that there is no UV divergence 
(due also to one dimension) . Therefore one can retrieve gauge invariant results by sending 
the UV cutoff to GO as is checked explicitly for the bosonic theory [17]. 

^^From the bulk point of view, once tlie near-horizon limit is cissumed, the only spatial length scale 
is AX ~ [qsNY^^Is, which corresponds on the boundary to the temporal scale At ~ {gsN)~^^^£s ~ 
(Pym^)~^^^- These estimates of AX and At are consistent with the space-time uncertainty relation 
AtAX>eS [25]. 
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Let us describe how we introduce the UV cutoff more in detail. First we take the 
static diagonal gauge A{t) = diag(ai, • • • ,aAr), where Ua can be chosen to satisfy the 
constraint maX(j(aa) — mina(aa) < 27r by using the large gauge transformation with a 
non-zero winding number. We have to add to the action a Faddeev-Popov term 



Spp = -^21n 

a<b 



sm ■ 



(3.2) 



and the integration measure for Oa is taken to be uniform. Having fixed the gauge com- 
pletely by the above procedure, we can introduce a cutoff A in the Fourier-mode expansion 



iumt 



iLunt 



(3.3) 



n=— A 



n=— A 



where a; = 27r//3, taking into account that we impose periodic boundary conditions on both 
Xf'^{t) and ip^it)- Using a shorthand notation 



fciH \-kp=n 

we can write the action (3.1) as 5" = + where 



(3.4) 



1 ^ 

n=— A 
■ A 

E 



(Xg — at 



.n=— A 



I rwj — 



(3.5) 
. (3.6) 



Integrating out the fermionic variables, one obtains the Pfaffian Pf A^, which is complex 
in general. (See eq. (B.19) for the definition of M.) The phase of the Pfaffian turns out to 
be quite small for SU(2) in the parameter region investigated in the present work, but it 
does fluctuate for SU(3), in particular at large /?. Here we simply neglect the phase and use 
\P{M\ = det(P^/''), where V = M'^M, instead of PfA^. The system with finite degrees 
of freedom we arrive at in this way can be simulated by using the Rational Hybrid Monte 
Carlo (RHMC) algorithm [36], which has become quite standard nowadays in lattice QCD 
simulations with dynamical quarks. The details of the algorithm are given in appendix B. 

An important assumption in our method is that the phase of the Pfaffian can be ne- 
glected. While the previous results [18, 19, 20, 24], which confirmed the gauge-gravity 
correspondence with high accuracy, certainly support this assumption, we cannot provide 
purely theoretical justification at the present stage. The effect of the phase can be incor- 
porated in principle by reweighting when one calculates expectation values. This direct 
method becomes impractical due to huge cancellations when the phase fluctuates violently. 
We suspect that the fluctuation of the phase is actually smaller than that of typical ex- 
tensive quantities at large /3 and that the phase-quenching is completely justifiable in the 
large-/? limit. In appendix C we investigate the effect of the phase on a typical observable 
and show that it is indeed negligible. 
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3.2 results for supergravity modes 

In this subsection we present our results for the two-point correlation functions of operators 
(2.18)-(2.20) corresponding to supergravity modes. 

Let us first define the operators we study by simulating the model (3.1). The operators 
Jtii-,ie (^>l)are defined by 

jtiL,H^^^^{^J^n---Xi^ ' (3-7) 

where F^j = —i [Xi,Xj] and Str represents the symmetrized trace treating F^j as a single 
unit. Restricting ourselves to £ <7, we may assume that all the indices i, j, ii, ■ ■ ■ , ie are 
different from each other so that the traceless condition for (2.18) is trivially satisfied. The 
operators T^j^_._ {£ > 2) are defined by 

rj:...,,-^Str(x,...X,,) . (3.8) 

Restricting ourselves to ^ < 9, we may assume that all the indices i, j, ii, i£ are 
different from each other so that the traceless condition for (2.20) is trivially satisfied. The 
operator T^^^ are defined by 

^il,..,, ^ ^Str((AX.)X., •••X,,) , (3.9) 

where Str represents the symmetrized trace treating (DfXi) as a single unit. Restricting 
ourselves to £ < 8, we may assume that all the indices i, ii, are different from 

each other so that the traceless condition for (2.19) is trivially satisfied. Due to the S0(9) 
symmetry, the result for the two-point correlation function for each type of operators 
should not depend on the assignment of the indices. Therefore we average over all possible 
assignment to increase the statistics in actual calculation. 

Since the basic dynamical degrees of freedom in our Monte Carlo calculations are 
Fourier modes (3.3), the correlation functions that arc directly accessible are those in the 
momentum space (^0(j)) 0(—p)'^, which are related to those in the real space by the inverse 
Fourier transformation 

{0{t) O{0)) = I ^ (dip) 6{-p)) e^^' . (3.10) 

Using the Fourier modes On defined similarly to (3.3), we can rewrite it as 

(0(p) 6{-p)) = 13 (On O-n) , where P = ^ • (3-11) 

Note that the factor of (3 on the right-hand side is needed to make the correlation function 
finite in the large-^ limit. 

The gauge-gravity correspondence predicts the two-point function in the Fourier space 
to behave as 

(0{p) d{-p)) ~ f{p) + g{p) \pf'' (3.12) 
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at small p as one can see from eqs. (2.14) and (2.15). Here, f{p) and g{p) {g{0) / 0) 
are analytic functions invariant under p -H- —p, and hence they can be written as f(p) = 
/o + /2 + ■ ■ ■ and g{p) = go + g2P^ + ■ ■ ■ , where 50 7^ 0. The coefficients are not fixed 
since we have focused on the most relevant term in the IR limit neglecting the overall 
normalization factor. Our main task is to extract the power u by fitting the Monte Carlo 
data for various correlation functions to the form (3.12) and to compare it with the values 
predicted by the gauge-gravity correspondence. 

The existence of the undetermined analytic terms in (3.12) makes the extraction of 
u more difficult, in particular, for large u. This problem can be avoided if one can make 
the inverse Fourier transformation (3.10) numerically. The analytic terms are transformed 
into local terms, and hence do not affect the power-law behavior (2.16) of correlation 
functions in the real space, from which one can extract the exponent u. However, there 
arc cases in which the inverse Fourier transformation is not possible numerically. This can 
happen either due to the IR behavior or due to the UV behavior. Correlation functions 
for operators with u < [e.g., T2^) behave as ~ \p\'^^ at small p, and hence the inverse 
Fourier transform is divergent at p ~ 0. The correlation function for operators {e.g., T^) 
including a derivative does not fall off in the momentum space at large which makes 
the inverse Fourier transformation numerically unstable. Even in these cases, it turns out 
that the results we obtain directly in the momentum space are in good agreement with the 
prediction from the gauge-gravity correspondence. 

In what follows we present our results for J^, T^^ and in order. In the N = 2 
case the correlation function becomes identically zero for J+ (£ :even), r++ (£ :odd) and 

{i : even) due to properties of the Pauli matrices, and hence we omit these cases. 

Let us start with the correlation functions of operators , for which the inverse 
Fourier transform can be calculated numerically. If we naively make the inverse transform, 
however, the correlation function in the real space shows oscillating behavior with the 
period St ~ /3/(27rA). This is well-known as the Gibbs phenomenon, and it is an artifact 
of the sharp cutoff in the momentum space. In the present case, we may naturally expect 
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Figure 1: The log-log plot of the correlator (^J/(t) ^7(0)^ with ^ = 1,3 for = 2 (Left) and 
with £ = 1, 2, 3, 4 for N = 3 (Right). The cutoff parameters are chosen as /3 = 4 and A = 16. The 
straight lines represent the power-law behavior predicted by the gauge-gravity correspondence. 
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Figure 2: The correlator (^Ji{p) Jii—p)"^ is plotted as a function of p for = 2 (Left) and = 3 
(Right). The extrapolation to A = oo is made at each p. The dashed line represents a fit to the 
behavior (3.15), from which we obtain u = 0.42(6) for A/' = 2 and u = 0.43(5) for N = 3. The 
predicted value is = 2/5. 



that the large-p behavior of the two-point correlation functions is given by 



J/(p)JH-p)}-^- (3-13) 

Indeed our data can be nicely fitted to this behavior at large p, and we can obtain the 
coefficient k reliably. The correlation function in the real space can then be obtained as 

= (j/o4) +E2COS (^) (j/J^-n) . (3-14) 
n>o \ y / 

where we extend the sum over n beyond A using the form (3.13) up to n = 1000. The 
results obtained in this way are shown in fig. 1 for A?^ = 2 (Left) and N = 3 (Right). The 
cutoff parameters are /3 = 4 and A = 16. We fit our data points for n = 10, 11 and 12 to 
(3.13), and use the form for 12 < n < 1000 in the sum (3.14). Straight lines represent the 
power-law behavior const. /t^^^"*^ with u = 2£/5, which is predicted for N = oo from the 
gauge-gravity correspondence. In the large-t region {t > 0.5), the power law is reproduced 
with remarkable precision up to the IR cutoff scale {t < /3/2). Moreover, the results for 
N = 2, which are obtained only for odd £ for the reason mentioned above, already show 
the power-law behavior with the same exponent as predicted by the supergravity analysis. 
This suggests that the exponents are actually independent of A'"(> 2). 

The consistency with the gauge-gravity correspondence can also be seen directly in 
the momentum space. Since we are interested in the IR behavior, let us restrict ourselves 
to the behavior at small p. According to the gauge-gravity correspondence, the first few 
terms relevant at small p for ^ = 1, 2, 3, 4 are 

~T+r^J+r A /a + 6H2- + cp2 for£=l,2, 
4iP) 4i-P)) ^[a + bp^ + a Ipl- for ^ = 3, 4 , ^'-''^ 

where u = 2i/5. In fig. 2 we show our results for with N = 2 (Left) and A^ = 3 (Right). 
At each p we made an extrapolation to A = oo assuming that the finite A effect is of the 
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h 


c 
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'^pred 




0.40(1) X 10-2 


-0.8(1) X 10-3 


0.32(8) X 10-4 


0.43(5) 


0.4 


4 


0.438(1) X 10-3 


-0.58(7) X 10-4 


0.26(9) X 10-4 


0.84(3) 


0.8 


4 


0.47 X 10-4 


-0.38 X 10-5 


0.17 X 10-5 


1.18 


1.2 


4 


0.582(6) X 10-5 


-0.12(5) X 10-6 


0.57(4) X 10-^ 


* 


1.6 



Table 2: The parameters obtained by fitting our results for the correlators (^J^{p) J^{—p)j with 

iV = 3 to the behavior (3.15). The values of v predicted from the gauge-gravity correspondence are 
shown on the right most column. For Jg", the fitting errors are not shown since we can use only 
4 data points for the fit. For , wc present the values of a. h and c obtained by fixing u to the 
predicted value 1.6 since the fitting errors turned out to be large otherwise. 
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Figure 3: The log-log plot of the correlator {f++{p)f++{-p)j for TV = 2 (Left) and = 3 
(Right). The extrapolation to A = oo is made at each p. The dashed line is a fit to the behavior 
(3.16), from which we obtain v = —0.57(3) and u = —0.66(2) for A'' = 2 and N = 3, respectively. 
For the N = 3 case, we can fit the data up to p/ (27r) = 2 as shown by the solid line by including a 
higher order term ^1^^"+^ in (3.16), which gives v = —0.61(2). The predicted value is v = —3/5. 



order of 0(1/ A).^^ By fitting our results at < p/(27r) < 1 to the behavior (3.15), we 
obtain v = 0.42(6) and v = 0.43(5) for iV = 2, 3, respectively, which are consistent with the 
predicted value u = 2/5. The same analysis can be performed also for with £ = 2,3, 4. 
The values of the parameters obtained by fitting our data for iV = 3 are summarized in 
Table 2. The fits are reasonable, and they are consistent with the predictions from the 
gauge-gravity correspondence. 

Let us next consider the operator for ^ = 2,3, 4, 5. The gauge-gravity correspon- 
dence predicts v = {2£ — 7)/5 [11]. The ^ = 2, 3 cases are of particular interest since the 
supergravity analysis predicts v = —3/5 and v = —1/5, respectively, which are negative. 
Prom (3.12), the correlation function in the momentum space is expected to behave as 

(T++ip) f ++(-p)^ = a IpI^" + b for ^ = 2, 3 (3.16) 

^®In most cases we used the data for A = 8, 12, 16. In the N = 3 case with p/(27r) < 3/0 ioT > 5, 

we had to use the data for A = 6, 8, 12 since the data for A = 16 had too large statistical errors. The 
extrapolation to A = (X) seems to be fine, though, since finite A effects are less severe in the small p region. 
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Figure 4: The log-log plot of the correlator (t^~^{p) f^~^{—p)'^ for N = 3. The extrapolation to 

A = oo is made at each p. The dashed line represents a fit to the behavior (3.16) within the range 
is 2/13 < p/{2it) < 0.8, which gives u = -0.2(1). The predicted value is u = -1/5. 
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Figure 5: (Left) The correlator ^T/+(p) T/+(-p)^ is plotted as a function of p for ^ = 4 (Left) 
and 1 = 5 (Right) in the N = 3 case. The extrapolation to A = oo is made at each p. The dashed 
line represents a fit to the behavior (3.17), which gives f = 0.19(3) and v = 0.7(1) for £ = 4 and 
£ = 5, respectively. For the £ = 5 case, we can fit the data up to p/{2it) = 2 as shown by the solid 
line by including a higher order term in (3.17), which gives v = 0.60(9). The predicted 

values are v = 1/5 and v = 3/5 for £ = 4 and 1 = 5, respectively. 



at small p. Note that the leading term is divergent in the IR limit p — t- unlike the previous 
cases. In particular, the inverse Fourier transform (3.10) is ill-defined for i = 2. This is 
somewhat reminiscent of the free-field behavior ~ 1/p^ albeit with a milder form. In fig. 3 
we show a log- log plot in the case of ^ = 2 for AT = 2 (Left) and N = 3 (Right). For 
the N = 2 case, we can fit all the data in the figure with p/(27r) > 1/P to (3.16), which 
gives u = —0.57(3). For the A'" = 3 case, we can fit the data within 1//3 < p/{2it) < 1 
as represented by the dashed line, which gives v = —0.66(2). By including a higher order 
term (i|pp'^+^ in (3.16), however, we can fit all the data in the figure with p/(27r) > 1/(3 
as represented by the solid line, which gives v = —0.61(2). These results are consistent 
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'^pred 




0.249(6) X 10-1 


-0.13(2) X 10-2 
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-0.61(2) 


-0.6 


rn++ 

-'s 


0.23(5) X 10-2 


-0.8(6) X 10-3 
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-0.2(1) 


-0.2 


rp+ + 

-'4 


0.38(3) X 10-3 


-0.16(3) X 10-3 


0.45(7) X 10-6 


0.19(3) 


0.2 


^5 


0.12(5) X 10-^ 


-0.17(3) X 10-5 


0.2(1) X 10-6 


0.60(9) 


0.6 



Table 3: The parameters obtained by fitting our results for the correlators {T^^ {p)T^'^ {—p)j 
with = 3 to (3.16) for = 2,3 and to (3.17) for £ = 4, 5. The values of v predicted from the 
gauge- gravity correspondence are shown on the right most column. For i = 2 and 5, we also include 
a term with the coefficient d = 0.95(6) x IQ-^ and d = -0.24(7) x IQ-^, respectively. 



with the predicted value v = —3/5. The deviation from the behavior (3.16) at smaller p 
around the IR bound (~ /3~^) is naturally interpreted as finite /3 effects. In fig. 4 we plot 
the result for Tg""*" in the N = 3 case. Fitting our data to the behavior (3.16), we obtain 
u = —0.2(1), which is consistent with the predicted value u = —1/5. 
In the ^ = 4, 5 cases, the gauge-gravity correspondence predicts 

(r++(^,)r++(-p)) =a + 6|^.|2- + cp2 (3.17) 

at small p, where v = 1/5,3/5 for i = 4,5, respectively. In fig. 5 we plot the correlation 
functions for T^~^ (Left) and Tc^'^ (Right) with = 3. The data at small p have large 
statistical errors. We consider this as a consequence of the fact that the operators T^"*"^ 
with large £ are composed of large powers of X without derivatives or commutators, and 
hence they are more affected by the large fluctuations of the low-momentum modes of the 
X field. For the N = 2 case, we can fit all the data in the figure to (3.17), which gives 
u = 0.19(3). For the = 3 case, we can fit the data with p/(2tt) < 1.2 as represented 
by the dashed line, which gives = 0.7(1). By including a higher order term djpp'^'''^ in 
(3.17), however, we can fit all the data in the figure as represented by the solid line, which 
gives u = 0.60(9). These results are consistent with the predicted values u = 1/5,3/5 for 
^ = 4, 5, respectively. In Table 3 we present the values of the fitting parameters obtained 
for r++ {e = 2, 3, 4, 5) in the AT = 3 case. 

As the last example of operators corresponding to supergravity modes, let us consider 
the operator defined by (3.9). The peculiarity of this operator is that it involves a 
derivative with respect to t. As a result, the corresponding two-point correlation functions 
do not decrease as ~ 1/p^ unlike the previous examples. Although finite A effects turned 
out to be rather large, we were able to extrapolate our results obtained at A = 6, 8, 12, 16 
to A = DO. In fig. 6 we plot (T}{p)f}{-p)^ for N = 3, 13 = 5. Fitting the data in the 
region < p/(27r) < 1 to the behavior 

(^f+{p)f^{-p)^=a + b\p\"' + cp' , (3.18) 

we obtain ly = 0.80(3), which is consistent with the predicted value u = 4/5. 

Finally, let us briefly comment on the sign of the leading term in all the correlation 
functions studied in this paper. Our data together with the formula (2.17) shows that the 
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Figure 6: The correlator (^T^ (p) {—p)^ is plotted for A'' = 3 and ^ = 5. The extrapolation 

to A = oo is made at each p. The dashed line represents a fit to the behavior (3.18) with a — 
2.614(2) X 10--\ b = 3.1(3) x 10-^ c = -1.1(3) x IQ-^, and u = 0.80(3). The value of u predicted 
from the gauge-gravity correspondence is v = 4/5. 

correlation functions obey the positive IR behavior in the real space suggested from the 
(reflection) positivity, except for the operators T^'^,T^ with £ = 2. In the case of T2"^, 
we note that the naive integrability condition for the inverse Fourier transformation is 
violated in the IR region since (T^'^{p)f++{-p)'^ ~ at p -> 0. The non-positivity 

may be attributed to this property in analogy with the apparent violation of positivity that 
occurs for the free-field propagator (^X{t)X{0)'^ = —\t\/2 + const., which is related with 
the divergent momentum-space behavior l/p^ at p = 0. In the case of T^, the momentum- 
space correlation function does not decrease for p ^ oo as exhibited in fig. 6, and hence 
the UV behavior does not satisfy the integrability, either. Note, however, that the non- 
integrability at the UV region does not directly affect the leading IR behavior. Therefore 
the situation is rather unclear. 

3.3 results for stringy excited modes 

In this subsection we present our Monte Carlo results for operators corresponding to stringy 
excited modes. Let us first consider the BMN-type operator (2.25). The result (2.43) based 
the gauge-gravity correspondence suggests that the corresponding two-point functions be- 
have as 

(o4„(i) OiliO)) ~ aexp [-bt'/') , (3.19) 

up to a power behaved correction factor. This prediction was obtained in the plane-wave 
limit, which is justified when the angular momentum J and the wave number n are both 
large. 

In what follows we calculate the correlation function on the gauge theory side and 
compare the results with the prediction (3.19). For that purpose we make an inverse 
Fourier transformation as we did for the operator in the previous section. In order to 
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Figure 7: The two-point correlation function of the BMN-type operator (2.25) is plotted in the 
real space for J = 2 (Left) and J = 3 (Right) with n = 1. The matrix size is A'' = 3 and the cutoff 
parameters are 1/ (3 — 0.15 and A = 16. The curves represent the fits to the form aexp (— &t^) with 
c = 0.3, 0.6, 0.9. Fitting the data in 0.2<t< 1.5, we obtain c = 0.63(1) and c = 0.66(2) for J = 2 
and J = 3, respectively. 



avoid the Gibbs phenomenon associated with the sharp UV cutoff, we extrapolate our data 
to larger momenta assuming the form ^I'^^i ■ We take the UV cutoff to be A = 16, and 
determine the parameters k and by fitting the data at p = 27r/c//3 for k = 8, • • • ,12. 
Similarly to the case of T^~^ with £ = 4,5 shown in fig. 5, the two-point functions of 
the BMN-type operators for J > 4 have large statistical errors in the small p region, 
which propagates to the real-space correlator through the inverse Fourier transformation. 
Therefore, we restrict ourselves to J = 2, 3. 

Figure 7 shows the results for BMN-type operators with J = 2 and J = 3, respectively. 
The wave number n in the operator (2.25) is set to n = 1. Fitting the data within 
0.2 <t < 1.5 to the form aexp {-bt"), we obtain c = 0.63(1) and c = 0.66(2) for J = 2 
and J = 3, respectively, which are close to the value 3/5 predicted for large J and n. The 
deviations observed at large t (> 2) can be attributed to finite-/3 effects. The coefficient 
b is predicted to be proportional to n/J for large n and J with n :$> J. The values of 
b obtained from the fit by fixing c to be 3/5 are b = 4.70 and b = 5.40 for J = 2 and 
J = 3, respectively. This suggests that the coefficient b has more nontrivial dependence 
than 6 oc n/J for small n and J. 

Next we consider a non-BPS type operator of the form 

Oki = Tr(^FkiJ2{X,f^ , (3.20) 

where F^i = —i [Xk,Xi]. Note that this operator does not belong to the type of operators 
in (2.18) since the traceless condition is not satisfied. In fig. 8 we show our results 
for the two-point correlation function (^Oki{t) O^j^i{0)^ , which look surprisingly similar to 
the results for the BMN-type operators. The log-log plot shown on the right reveals clear 
deviation from a straight line behavior corresponding to the power law. Fitting the data in 
the region 0.2 < i < 1.5 to the form aexp {—bf^), we obtain c = 0.58(2), which is consistent 
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Figure 8: (Left) The two-point correlation function of the non-BPS operator Oki is plotted for 
N = 3, = 0.15, A = 16. The curves represent the fits to the form aex.p{—bt^) with c = 
0.3,0.6,0.9. Fitting the data in the region 0.2 < t < 1.5 gives c = 0.58(2). (Right) The log-log 
plot of the same correlation function, which reveals clear deviation from a straight line behavior 
corresponding to the power law. The solid line represents a fit to the behavior aexp(— 6/:^) with 
c = 0.6. 



with c = 3/5. This suggests the existence of a universal mechanism for generating a mass 
scale in the case of non-supergravity operators. 

4. Summary and discussions 

In this paper we have studied the strong coupling dynamics of Matrix theory by investi- 
gating the IR behavior of two-point correlation functions. As we have reviewed in section 
2, various predictions can be obtained from the bulk supergravity analysis based on the 
gauge-gravity correspondence. The direct calculations on the gauge theory side are com- 
pletely out of reach of perturbative analyses, however. This motivated us to perform 
detailed Monte Carlo evaluation of the gauge-theory correlators. Let us briefly summarize 
the salient features of our results. 

(a) The correlation functions of operators corresponding to the supergravity modes show 
the power-law behavior with the power consistent with the prediction. Note, in 
particular, that our Monte Carlo results are obtained at = 2 and 3, while the 
supergravity analysis is valid naively at iV = oo. Finite N effects, if they exist, would 
be of the order of which is 25% and 11% for AT = 2 and N = 3, respectively. 
Since the agreement with the supergravity predictions is much more accurate, it is 
conceivable that the power is actually independent of N {>2). 

(b) The power-law behavior seems to extend beyond the range (1.2) suggested by the 
validity of the supergravity analysis. Note that our results for AT = 2 shown in fig. 1 
(Left), for instance, include the data around t ~ 1.5, while the upper end of the range 
(1.2) is around t ~ 1.4 and 1.7 for N = 2 and 3, respectively. It seems reasonable 
to consider that the power law would continue to be valid beyond this bound if we 
increase the IR cutoff /3. Indeed, the results in the momentum space shown in fig. 2 
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suggest that the power law is vahd even at p/{2tt) ~ 0.2, which is much smaUer than 
1/t ~ (1.5)^^ ~ 0.72. This points to the interesting possibihty that the power law 
with the particular power predicted within the range (1.2) is actually valid even in the 
far IR region f ~ iV — > oo corresponding to M-theory. We consider it important to 
explore the implication of this statement on the wave function of the 11-dimensional 
graviton. In particular, the meaning of the anomalous behavior [14] with respect to 
the IID boost transformation discussed in section 2.3 must be clarified. 

(c) As for the stringy excitations, we only have semi-quantitative predictions for the 
BMN-type operators with large J and large wave number n. Our Monte Carlo 
results for the BMN-type operators with J = 2, 3 provide some evidence that the 
predicted form exp(— is valid even at small J, although the coefficient b has 
J-dependence different from what is expected at large J. The same form is also found 
to fit well our results for other non-BPS operators corresponding to stringy excited 
modes. Our results therefore suggest that this peculiar behavior may be a universal 
feature of the mechanism for mass scale generation in Matrix theory. 

In order to understand the physical meaning of the observed behaviors of the correla- 
tion functions, we have discussed the spectral representation of the correlation functions 
in appendix A. The power-law behavior for the supergravity modes can be naturally un- 
derstood as a consequence of the existence of the zero-energy bound state. On the other 
hand, the exponential behavior exp(— 6t^/^) for the stringy excited modes corresponds to 
a very peculiar form of the spectral function. It would be interesting to clarify possible 
implications of these properties on M-theory. 

Let us list some possible directions for future works. First of all, it is obviously worth- 
while to extend the Monte Carlo studies to larger N and larger /3, which would further 
clarify our observations (a) and (b) given above. It would be also important to study 
various other operators corresponding to stringy excited modes and to test the universal 
behavior stated in (c). On the technical side, it is important to clarify the issue related to 
the sign problem in Monte Carlo simulation discussed in appendix C. 

From the viewpoint of testing the gauge-gravity correspondence, it is also interesting 
to pay attention to the normalization, including the sign, of the correlation functions. 
For instance, the Ward-like identities for correlation functions are discussed [28] using the 
gauge-gravity correspondence in the non-conformal case including the case of DO-branes. 
It would be interesting to test such identities by our Monte Carlo method, which would 
also serve as a self-consistency check of the method itself. 

From a more general perspective of the gauge-gravity correspondence, it would be 
interesting to extend our work to gauge theories in higher dimensions. Of particular interest 
is to test the AdS-CFT correspondence by studying four-dimensional A/" = 4 SU(Ar) SYM 
by Monte Carlo methods, since most of the tests given so far are not quantitative except for 
some special cases in which perturbative analyses on the gauge theory side turned out to be 
useful. In fact, a first step in this direction is taken already. In refs. [37], matrix quantum 
mechanics of 6 bosonic commuting matrices [38] is shown to give results consistent with 
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the AdS-CFT correspondence for the three-point functions of chiral primary operators. 
More recently, Monte Carlo studies of the 4d = 4 SYM have been performed [39, 40] 
based on the novel large- A?" reduction [41], which reduces the calculation in the SYM (in 
the planar large- AT limit) to that in the BMN matrix model [29]. The advantage of this 
method is that supersymmetry is maximally respected. In particular, there is no need 
for fine-tuning the parameters in the regularized theory unlike the proposals based on the 
lattice regularization [42]. See ref. [43] for a proposal to combine a two-dimensional lattice 
with a fuzzy sphere in order to study the 4d theory at finite A^ without fine-tuning. 

It is also interesting to study the two-dimensional M = 8 SU(A) SYM by Monte Carlo 
methods. This corresponds to a system of coincident Dl-branes, which has an interesting 
connection [44] to the black hole/black string transition as exhibited by the phase diagram 
with respect to the temperature and the spatial volume. This system has been studied 
in the strongly coupled (low temperature) regime by using the dual type II supergravity, 
and the existence of a first-order transition is conjectured [44]. On the other hand, the 
weak coupling (high temperature) regime can be described effectively by a bosonic matrix 
quantum mechanics, in which two phase transitions of second and third orders are found 
[45, 46]. It is interesting to investigate how the two transitions merge into one as the 
coupling constant is increased (i.e., as the temperature is decreased) by direct Monte Carlo 
methods. A lattice simulation of two-dimensional SYM with 4 supercharges and the SU(2) 
gauge group was started a few years ago [47], and was extended to SU(A) with A = 2, 3, 4, 5 
[48]. It was shown that fine-tuning of parameters is not necessary in these cases. Recently 
the 16 supercharge case [49] has been studied. Let us also recall that the (1-|-1)D SYM 
with 16 supercharges is termed "Matrix String Theory" as a possible non-perturbative 
description of type IIA superstring theory in ten dimensions [50]. This theory is closely 
related with the doubly-compactified super membranes in 11 dimensions [51], and hence 
with Matrix theory studied in this paper. It would be very interesting to explore various 
non-perturbative relations between the two gauge theories by Monte Carlo methods. 

To conclude, with the various new ideas on how to treat supersymmetry on a computer, 
we consider that interesting non-perturbative physics in supersymmetric gauge theories, in- 
cluding those relevant to string theory, has become accessible by Monte Carlo simulation. 
We hope that our work not only provides further motivation to pursue the Matrix theory 
conjecture, but also encourages further development in Monte Carlo studies of supersym- 
metric gauge theories. 
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A. Spectral representation of the correlation functions 



In this section we discuss the physical meaning of the two-point correlation functions in 
Matrix theory suggested from the gauge-gravity duality and from direct Monte Carlo cal- 
culation. For that purpose we consider the spectral representation 



o(i)O(0)) = p d^ip{^i)l^{t,^i) , (A.i) 

in the Euclidean space. 

Let us first discuss the form (2.16) for the supergravity operators. This corresponds 
to the spectral representation p{n) ~ /Lt^'^+^. For 2i/ + 1 > 0, the //-integral in (A.I) is 

convergent and one retrieves (2.16). For 2z/ + 1 < 0, we again invoke analytic continuation 
in z^. This is in sharp contrast to the free-theory behavior (n = 1,2,3,- ••), which 
corresponds the spectral functions with delta-function-like singularities at = 0. The 
milder power-law singularity in the spectral function may be regarded as a signature for 
the existence of a zero-energy bound state as a result of very non-trivial IR dynamics. 
The fact that v increases with the angular momcntTim £ is consistent qualitatively with 
the existence of many-body states which consist of the zero-energy bound states, since the 
distance scale of such states can grow faster for larger i with a power-law behavior as the 
energy increases. 

Let us next discuss the operators corresponding to stringy excited modes. According 
to the scaling argument in sections 2.3 and 2.4, the general form of two-point correlators 
is given by 

0{t) 0(0)) ~ ,(^+«)/^ |t|-(7AH-12)/5 f . (A.3) 

For the BMN-type operators with large angular momentum J and a large wave number n, 
one obtains the prediction for the function /(5s-/V|ip/^^) at large distances as [16]^^ 



J 



apart from a possible power-law prefactor, with K being an operator-dependent constant. 
As the spectral function corresponding to the behavior (A.4) , let us consider 

p(/x)~exp(C/x«+i) , (A.5) 

up to a possible power-law correction factor. At large \t\, the //-integral is dominated by 

/ \ l/a 

the saddle point /t = ( |t|/(a -|- l)Cl . Comparing the resulting integral with (A.4), we 
can determine the parameters a and C. Plugging them into (A.5), we get 

p(/x)^exp^-^(:^) gV2^-3/A . (A.6) 




'^The time coordinate used in the first paper of ref. [16] is rescaled from ours by t — ^ q^^^t. 
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This is quite different from a simple delta- function p(/i) ~ — M), which represents the 
standard asymptotic behavior of the propag ator e-*^l*l for a massive particle state. At low 
energy — > 0, the spectral function (A. 6) decreases exponentially, which suggests that the 
single-body zero-energy bound state does not directly contribute to the intermediate states. 
On the other hand, the expression (A.6) approaches a constant if we naively extrapolate 
it to larger fi. However, it is expected that the power-law prefactor in (A. 3) becomes 
important in this regime. Since the prefactor is essentially the same as in the case of 
operators corresponding to supergravity modes due to the scaling symmetry, we consider 
that the many-body states consisting of the zero-energy bound state play a dominant role 
at high energy. 



B. The algorithm for Monte Carlo simulation 

In this section we explain the details of the algorithm we use for Monte Carlo simulation. 
In order to treat fermions efHciently, it is crucial to use the idea of the Hybrid Monte Carlo 
(HMC) algorithm [52]. 

Let us first discuss the idea for the bosonic part (3.5) neglecting the fermionic part. 
We introduce the (fictitious) momentum variables 

A 

Hf(t)= ^ n«?e*'^"*, (B.l) 

n=-A 

which are conjugate to X"''(i). Due to the Hermiticity of n"''(t), we have the constraint 

nf-n = (n'n)* ■ (B.2) 

Similarly, we introduce the momentum variables Pa corresponding to a^. Then we define 
the fictitious Hamiltonian for the bosonic part as 

N A 

= 2 T.(p-f + 2 E n^,-n + 5b[x, «] . (B.3) 

a=l n=— A 

The Hamilton equation can be obtained as 



dr 5n«'' 

dfl^ _ dHy, _ dSb[X,a 



ntn (B.4) 

(B.5) 



dr dX9-^ dX^^ 

in in 

By using the explicit form of S'b [X, a] , we can calculate the second line explicitly as 

dm' 



-Np[(nco-haa-a,)yxf% 



dT "LV /?' 

+2{X,X,X,t\-AiX,X,X,t\ + 2{X,X,Xif_\ . (B.6) 
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Similarly, the Hamilton equations for the gauge variables are obtained as 

^ = = 2iV X; (nc. - -{aa - a,))xtxf^r, ■ (B.8) 

n=— A 

Updating the bosonic variables is performed by (i) refreshing the momentum variables pa 
and by Gaussian random numbers, which obey the distribution c^^^ and (ii) solv- 
ing the Hamilton equations for a fixed time interval r. In order to satisfy the detailed 
balance, one needs to adopt the so-called leap-frog discretization for r-evolution, and to 
perform the Metropolis accept/reject procedure after (ii) with the acceptance probability 
mm 

(l,e-^^b). For more details, see refs. [53, 54], in which the method is described for 
supersymmetric matrix models. 

Let us next introduce fermions. We expand the Fourier component of the fermion 
field as 



N 



2 



A=l 



in terms of U(A'^) generators t . (We will impose the traceless condition later.) Here we 
use a representation given by 

where a = 1, • • • , N'^, and we have defined 



{i^)ab = SaiA^jA , (B.IO) 



A = N{iA-l)+jA , (B.ll) 
A = N{jA - 1) + u , (B.12) 



with which we may write 



Tr{t''t^) = 5AB ■ (B.13) 

Let us also define 

QABC = Tr = Sj^igdjgi^dj^i^ - Sj^igdj^i^dj^i^ . (B.14) 
Then the fermionic action S{ given in eq. (3.6) may be written in the form 

where we have defined a matrix M' of dimension 16(2A -|- 1)N^ by 

■M'Aap,BISq = (7i)"^ QABC + ^ Spq (ipio - i^{aig - ajj^ d^f} Sab ■ (B.16) 
We impose the traceless condition on the fermion field by making the replacement 

N-l 

{^anf" ^ -Y.^^-nr . (B.17) 
a=l 
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The fermion action can be written in terms of the remaining degrees of freedom as 

^ap M'Aap,BPq ^fq = ^ap MA'ap,B'^q i^f'q , (B.18) 

where we have defined a matrix M. of dimension 16(2A + 1)[N'^ — 1) by 

MA'ap^B'Pq = M'A'ap,B'l3q " ■^'N^ap,BI3q^iA' 3 A' 

~-^'A'ap,N2(Sq^iB'jB' ~^ ■^'N'2ap,N'^l3q^iA'jA'^iB'jB' ' (B.19) 

Integrating oTit the fcrmions, we obtain PfA^, which is complex in general. As we explained 
below (3.6), we simply replace it by its absolute value 



det 



(P^/^) , (B.20) 



where V = M^M. 

The trick of the RHMC algorithm [36] is to represent \PiM\ as 



PfA^ 



dFdF*e-^^^ , (B.21) 



where 

Q 

SpF = aoF*F + ^akF*{V + bk)-^F , (B.22) 

fc=i 

using the auxiliary complex variables F, which arc called pscudo- fcrmions. Here the con- 
stants a/j and are real positive parameters appearing in the rational approximation 



X 

k=l 

and they can be generated by a code [55] based on the Remez algorithm. The approximation 
(B.23) can be made to have sufficiently small relative errors (smaller than 5) within a certain 
range (e < a; < 1). In our simulation, we use Q = 15 with S = 1.19 x 10~^ and e = 10"^^. 
We rescale the matrix T> by an appropriate constant factor so that the largest eigenvalue 
of V is well below the upper bound for the approximation. 

Then we apply the usual HMC algorithm to the whole system. The Hamiltonian for 
the pseudo-fermions is given by 

FpF = ^^>l^>/ + SpF , (B.24) 
I 

where the field $ represents the momentum variables conjugate to the pseudo-fermions F. 
Here and henceforth, the index / is used to represent the spinor index a, the momentum 
index n and the SU(A^) index A, collectively. The Hamilton equations for the pseudo- 
fermions are given by 

dr d^i 

^ k=l 
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where G\ ' is defined by 



(m^'M + bk)^^ of = Fi for A; = 1, • • • , Q 



(B.27) 



We have to add extra terms 




(B.28) 



(B.29) 



on the right-hand side of eqs. (B.5) and (B.8), respectively. 

The main part of the computation comes from solving a linear system (B.27). We 
solve the system for the smallest hk using the conjugate gradient method, which reduces 
the problem to the iterative multiplications of to a pseudo-fermion field, each of which 
requires 0{I\?N'^) arithmetic operations. The solutions for larger ^^'s can be obtained 
as by-products using the idea of the multi-mass Krylov solver [56]. This avoids the factor 
of Q increase of the computational effort. 

When we solve the fictitious classical Hamilton dynamics, the step size of the dis- 
cretized evolution may depend on the Fourier mode. We take the step size to be pro- 
portional to the average fluctuation of each mode of the bosonic matrices so that tlic 
configuration space can be swept out most efficiently. This technique is called the Fourier 
acceleration [35] . In the lattice gauge theory, in order to apply the Fourier acceleration one 
has to transform the configuration to the momentum representation. On the other hand, in 
the present momentum cutoff method, the Fourier acceleration can be implemented with- 
out any additional cost since we are working directly in the momentum space. Thanks to 
this advantage, the efficiency of the algorithm is enhanced drastically. 

In actual simulation, we observe certain instability [18], which is related to the existence 
of the fiat direction [Xj,Xj] ~ in the potential term of the action (3.1). This instability 
can be seen by probing the observable 



Typically this quantity fluctuates around some value, but sometimes it grows rapidly and 
the simulation bumps. In order to avoid this problem, we introduce a cutoff on B?, which 
is taken to be sufficiently larger than the upper edge of the fluctuation. The value of the 
cutoff is 5.0 and 3.8 for A'" = 2 and N = 3, respectively. 

C. Sign problem 

In this section we discuss the so-called sign problem in Monte Carlo studies of Matrix 

theory (3.1). As we explained below eq. (3.6), we neglect the phase of the Pfaffian that 

appears from integrating out fermionic matrices. 

'^^Note that one should not construct M expUcitly and multiply it to a pseudo-fermion field literally, which 
would require 0{A^^N'^) arithmetic operations. Instead one should reduce the procedure to multiplication 
of N X N matrices using the original definition of M. 




(B.30) 
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Figure 9: (Left) The VEV (cos2r)o,a; is plotted for SU(3), /3 = 6.67, A = 8. (Right) The 
distribution ^('^^(a;) and the product p{x) <x p^°\x)'w{x) with the normaUzation / dxp{x) = 1 are 
plotted by the solid line and the dashed line, respectively. SU(3), 13 = 6.67, A = 8. 



Here we present some numerical evidence, which shows that the effect of the phase is 
indeed small. As a typical observable, let us consider defined by (B.30), and discuss 
how it is affected by the fluctuation of the phase. For that purpose we use the factorization 
method [57, 58, 59], which was proposed to study the effect of the phase in a general model 
suffering from the sign problem. Let us define the distribution function 



p{x) = (^6ix - B?) 
p(o)(x) = {5{x-R^) 



(C.l) 
(C.2) 



for the full model and for the phase-quenched model, respectively. One can then easily 
show that 

p{x) p^^\x)w{x) , (C.3) 



where the correction factor w{x) is given by 



w{x) 



I 



0,a; 



cosT 



0,x 



(C.4) 



Here, T represents the phase of the Pfaffian, and the symbol ( • )o,x represents a VEV with 
respect to the phase-quenched model with the constraint = x. 

Since the calculation of PfA^ is time-consuming, we calculate dctA^ = (PfA^)^, from 
which we can obtain w{x) = ( cos(2r) ) . This is sufficient for estimating an upper bound 

\ / Q,x 

on the effect of the phase since the factor of 2 in the cosine only magnifies it. In fig. 9 (Left) 
we plot w{x) iov N = 3, /3 = 6.67 and A = 8. In fig. 9 (Right) we plot the distribution 
function (x) for the phase-quenched model and the product p{x) oc (x) w{x) with the 
normalization J dx p{x) = 1. The difference between p{x) and p^^\x) is indeed negligible, 
which implies that the expectation value {R'^) is not affected by the effect of the phase. 

We speculate that the phase quenching can be completely justified in the large-/3 limit. 
At large /3, the expectation value of (i?^) can be obtained by solving the saddle-point 
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equation 

-^lnp(o)(a;) = —^lnw(x) . (C.5) 
ax ax 

First, we can easily prove that the Pfaffian becomes real if we omit the kinetic term in 
(3.6). As P increases, there are actually more and more low-momentum modes, for which 
the kinetic term is small. Therefore, we consider it conceivable that the fluctuation of 
the phase docs not increase as 0{\/^), which is the typical growth of the fluctuation for 
extensive quantities at large (3. As a result, the right-hand side of (C.5) is expected to 
be 0(/3P) with p < I, whereas the left-hand side is expected to be 0(/3) according to the 
usual scaling argument. In fig. 9 we do observe that the x dependence of w{x) is much 
smaller than that of p^^\x). If our scenario is correct, we can neglect the effect of the phase 
completely in the /3 — t- oo limit. In order to confirm this scenario, we need to investigate 
the P dependence, which we leave for future investigations. 
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